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Abstract 



We calculate the total and differential decay rates for the semi-inclusive non- 
leptonic decays B — > D s X g (q = c or u) under the factorization hypothesis. 
The initial bound state effect is treated using the light-cone expansion and 

the heavy quark effective theory. We investigate the contribution of the pen- 

(*) 

guin amplitude for the decay mode B — > D s X c and find that it is small 

but non- negligible. The resulting decay rates for B — > D s X c agree with the 

(*) 

measurements. The use of the decay mode B — > D s X u to determine \ V u b\ is 
discussed. 
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I. INTRODUCTION 



The semi- inclusive nonleptonic decays B — > Dy>X q (q = c or it), where X g is a hadronic 
system containing a g-quark, can be used to study mechanisms of the production of the 
D(*> meson in B meson decays, as well as the dynamics of the strong interactions [|TJ. 
As pointed out in Ref. 0, the Cabibbo-Kobayashi-Maskawa (CKM) matrix element \V u b\ 
can be determined by the measurement of the momentum spectrum of the meson in 
B — > Dy'X u . An advantage of the method^ is that the majority of the spectrum in the 
rare decays B — > D^*>X U lie above the kinematic limit for B — > D^X C which constitutes an 
overwhelming background. Since this method involves nonleptonic B decays, in which only 
hadrons appear in the final state, the hadronic complications may be more sever than the 
determination of \ V u b\ from semileptonic decays B — > lviX u . An understanding of the strong 
interactions, which are responsible for the confinement of quarks and gluons in hadrons, in 
the weak decays is essential to a determination of \V u b\. 

Measurements of the branching fraction for B — ► D S X have been reported by ARGUS 
and CLEO. The average of their results gives B(B -> DfX) = (10.0 ± 2.5)% @. Recently, 
new precise measurements of and D* + meson production from B meson decays have 
been presented by BABAR ||. The results for the branching fractions are B(B — > D+X) = 
(10.93 ±0.19 ±0.58 ±2.73)% and B(B -> D*+X) = (7.9 ±0.8 ±0.7 ±2.0)%. As for the rare 
decays B -> D^X U , none of them have been experimentally measured. With high-statistics 
event sample from the B factories, coupled with improvements in detector performance and 
analysis techniques, measurements of B — > D^X U might be possible in the near future. 

The b — > qcs tree transition is expected to dominate B — > D^~X q decays. In these 
processes the b quark decays to a c or a it quark, emitting a virtual W that fragments to a 
cs quark pair that subsequently hadronizes into a meson. The effective weak Hamil- 

tonian responsible for the decays is well known ||, which is based on the operator product 
expansion [0 and the renormalization group ||. The problem is how to calculate hadronic 
matrix elements of four-quark operators in the effective weak Hamiltonian. Factorization 
hypothesis has widely been used in calculations of nonleptonic B decays, which reduce the 
hadronic matrix element of a four-quark operator to the product of two matrix elements of 
current operators 0. The argument for factorization is based on the space-time evolution 
of the decay products [|nj. At the quark level the decay B — > D^*'~X q begins as a nearly 
collinear configuration of b — > q ± (cs), with q and the (cs) pair moving rapidly apart in 
opposite directions. When the color singlet -D^ - meson is formed, it is far away from the 
hadronic system X q recoiling against it. The strong interaction between them is expected to 
be not significant. Therefore, the hadronic matrix element can be factorized into a product 
of hadronic matrix elements of color-singlet quark currents. For B — > D^~X C due to limited 
energy release this argument is considered to be weaker than B — > D(*>~X U . Factorization 
in the nonleptonic B decays with the emission particle being a light meson can be justified 
more rigorously in the heavy quark limit using perturbative QCD, where the emission par- 
ticle means the one that does not inherit the spectator quark from the B meson |TT|--|T3"|| . 



1 Other methods of extracting \ V u b\ from nonleptonic B decays have also been proposed ||. 
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Unfortunately, this is not the case for both B -> D^~X C and B -> DW-1 M decays, where 
the emission particle is that is not light. Nevertheless, factorization has been found 

to hold with current accuracy in the color-favored decay modes B° — > , where no 

perturbative QCD justification has been presented [Q. Thus, we have some confidence in 
the applicability of factorization to the color-favored decay modes B — > Z)W - X g . 

In this paper, we analyze 5 — > D^~X q decays in some detail. We shall calculate the 
total and differential B — > D^*'X q decay rates assuming factorization. We include the 
contributions of penguin diagrams to B — > Dy>~X c decays, which turn out to be small but 
non- negligible. After factorization the long-distance QCD effects on B — > Dy'~X q decays 
are contained in two matrix elements: One is related to the decay constant for the Dy^ 
meson, accounting for the direct generation of the Dy*> meson from the vacuum by the cs 
current forming the weak Hamiltonian, while the other incorporates the effect of the initial 
b quark bound state in the B meson, which relates a 6-quark decay to a 5-meson decay. 
The treatment of the initial bound state effect will be the main focus of this paper. We shall 
use the light-cone expansion and the heavy quark effective theory (HQET) jnj to make our 
calculations aimed at understanding this effect from QCD. 



II. FACTORIZATION 

Both B — ► Dy*^X c and B — > Dy*^~X u involve tree amplitudes. In addition, penguin 
amplitudes contribute to B — > Dy'X c . The annihilation and exchange diagrams are esti- 
mated to give much smaller contributions than the tree diagram || . We shall neglect them. 
The relevant AB = 1 effective Hamiltonian reads 

n eS = -| lv qb v^( Cl oi + c 2 oi) - v tb v* -£c n o n \ , (i) 

where Of 2 are the tree operators and 3 _ 6 (O 7 _i ) the QCD (electroweak) penguin opera- 
tors. They are given by 

0\ = (qibi) V -A(sjCj)v-A, 0\ = (qibj)v-A(sjCi)v-A, 

3 (5) = (Sibi)v-Aj2%<lj)v-(+)A, °4(6) = (Sibj) V -A J^i^'i) V-(+)A, 
q> q> 

3 3 
07(9) = 2& h i)v-AY^ e *Wj)v+(-)Ai °s(io) = -{sibj) V -Aj2e q >%ql)v+(-)A, (2) 

q' q> 

where % and j are color indices, (<?i<?2)v±a = 9i7^(l ± 75)<?2, the summation runs over 
q' = u, d, s, c, b, and e q i is the electric charge of the q' quark in units of \e\. 

The Wilson coefficients c n have been calculated in different schemes 0. In this paper 
we shall consistently use the naive dimensional regularization scheme. Including the next- 
to-leading order QCD corrections, the values of c n at the renormalization scale m b are 

c x = 1.082, c 2 = -0.185, c 3 = 0.014, c 4 = -0.035, c 5 = 0.009, c 6 = -0.041, 
c 7 = -0.002a em , c 8 = 0.054a em , c 9 = -1.292a em , c w = 0.263a em . 

Here a em = 1/137 is the electromagnetic fine structure constant. 
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Under the factorization hypothesis, we obtain from the effective Hamiltonian in Eq. (|]) 
the decay amplitudes 



A(B - D-X c ) = if Da (ap% s (X c \jC\B) + (3(X C \J\B 
A(B - D*-X c ) = am D *J D *e>**(X c \j;\B), 

i^v ub v; sai f D y Da (x u \j;\B), 



A(B -> D~X U ) 



G 



where 



A(B -> = -|K 6 K>im D */D*e^(X u |j;|S), 



a = ^[KbV^ai - K 6 y t *(a 4 + ai ) 



/3 = ^^:(a 6 + 



2m 2 



«2i 



C2i + 



C2i-1 

IvT 



tt2i-l — c 2i-l + 



C2i 



(3) 
(4) 

(5) 
(6) 



(7) 

(8) 
(9) 
(10) 



N c is the number of colors. The currents = 97^(1 — 7s)& and J = c(l — 75)6. 

We have defined the decay constants of the pseudoscalar D s meson and the vector D* s 
meson as 



(D7(p Ds )\sYrc\Q) = ifnA 
( D t~(Prn» e )\^ c \ Q ) = m D *f D *e 



(A)* 



(11) 
(12) 



where e stands for the polarization vector of the D* meson. Using the definition ([11]) and 
the Dirac equation, it follows that 



(£>7|s(l + 75 )c|0)=i/ 



' m c + m s 



(13) 



which has been used to obtain Eq. @. 



III. INITIAL BOUND STATE EFFECT 

In this section we study the initial 6-quark bound state effect on B — > D^~X q using the 
light-cone expansion and the heavy quark effective theory. In the following we work out, in 
detail, the formulation for B — > D~X C . The results for B — > D* S ~X C and B — > D^~X U , 
which have a simpler structure in decay amplitudes, can be easily obtained along the same 
lines. 

The differential decay rate for B — ► DjX c in the B rest frame is given by 

dT{B -> D-X c ) = 1 f*°- J2(27r)W(p B - p Ds - p Xc )\A(B - D-X c )\\ (14) 
2m B {2Tiy i 2E Da 
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Applying Jd 4 y exp[-iy ■ (p B - p Ds - p Xc )] = (2tt) 4 5 4 (p b - Pd 3 ~ Px c ), (X c \Q^{0)\B) = 
(X c \Q^(y)\B)exp[—iy ■ (px c — Pb)} due to translation invariance, and the completeness of 
the hadronic states \X C ), we obtain from Eq. ([3]) 

£ (2n) 4 5 4 (p B -p Ds -p Xc )\A(B -> D-X c )\ 2 

x c 



= fk J dSe^{\a\Y D y D M[jt\^fM\B) + \P\ 2 (B\[j\0), J(y)]\B) 

+aPY D AmJ\o),fM\B) + ^Pp" D Amj^oiAy)]\B)). (is) 

Computing the commutators of currents yields 
E W5\p B -p Ds -p Xc )\A(B - D-X c )\ 2 

-2f Ds (\a\ 2 S, a ^ D y DB + \f3\ 2 g a ?) [ tfye»*>» [d a A c (y)] (B\b(0)^U(0,y)Ky)\B) 



x. 



-2f Ds (a(3* + a*(3) m y Ds J d'ye^zA c (y)(B\b(0)lpU(0,y)b(y)\B}, (16) 

where S mufi = g^ a g u p + 9np9va - 9^9ap and 

A q (y) = -^JL_ / d % e- ik -ve(k°)5(k* - mj). (17) 

The soft gluon interactions on the final state c quark appear in the Wilson line operator 

U(x,y) = Vexp[ig s [* dzfA^z)], (18) 
Jy 

where V denotes path-ordering and is the background gluon field. The initial 
bound state effect is contained in the matrix element of the non-local 6-quark operator, 
(B\b{G)7 p U{Q,y)b(y)\B). We note that the same matrix element also incorporates bound 
state effects in other inclusive £>-meson decays, including B — > ti>tX q |L6|], B — > ^X s |L7 



and B -> K&X Pfl. We shall use the light-cone expansion to calculate the matrix element. 
The light-cone expansion enables a systematic ordering of nonperturbative QCD effects as 
an expansion in powers of a small parameter of order Aq CD /m^. The general method has 
been described in |16|JT7|| . We shall restrict ourselves here to the essential steps only. 

Each integrand of the integrals in Eq. ( |T5| ) contains an oscillating factor e iy ' PDs . The 
dominant contribution comes from domains with less rapid oscillations, i.e., \y ■ Pd s \ ~ !■ 
This implies that for sufficiently large momentum of the outgoing D s meson, the decay 
dynamics is dominated by spacetime separations in the neighborhood of the light cone 
y 2 = 0. The light-cone dominance is also justified since the function A c (y) in the integrals 
in Eq. ( TIB] ) has a singularity at y 2 = 0. Therefore, at leading twist we have[| 



Higher twist effects have been analyzed quantitatively in [19;. They give rise to corrections of 
order Aq CD /m| and can be added systematically. The calculation of higher twist effects on the 
processes under consideration is beyond the scope of this paper. 
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{B\b(0)^U(0,y)b(y)\B) = 2p B f <K< '^/(O- (19) 

Jo 

where /(£) is the 6-quark distribution function of the B meson introduced in Refs. [16,|l7l 

m = i_j d(y^ e ^ fl{5 |5 (0) ^ (0)y)6(y) | 5) |^^ (20) 



Note that /(£) depends only on the properties of the i?-meson bound state and is not specific 
to particular 5-meson decay process in question. Thus, /(£) is a universal distribution as 
fundamental as parton distributions in deep inelastic scattering. It has the interpretation of 
the probability of finding a 6-quark with momentum £p B inside the B meson with momentum 

/'/; OH- 

Assembling all the pieces, we obtain the momentum spectrum for D s in B — > D s X c 

dT , 5 fk\pvf r 1 



IT En. JO 



'I !>/> " ' ' * E D S JO 

x5(m 2 B £ 2 - 2m B E Da £ + m 2 Dg - m 2 c ) 

x [\a\ 2 (m 2 Da E Da + m B m 2 D £ - 2m B E 2 D ^) + \(3\ 2 {E Da - m B £) 
-{a(3* + ct(3)m c E Ds }. (21) 

Carrying out the 5-function integration in Eq. fl2~ll) , we finally arrive at 

dT -(b^d;x^ fh ° IpdJ2 



d \PD s \ 2 * m BE Da y/\l> Da \ 2 + m? c 

x{f(e + )[\a\ 2 (2m B E 2 Ds e + - m B m 2 Ds t% - m 2 D E Ds ) 

+|/3| 2 (m B e; - E Ds ) + (a/3* + a*/3)m c E De ] 

-(£-£)}, (22) 



where 



= Ds v ' PDsl — (23) 

The interference between the tree and penguin amplitudes is associated with the coefficients 
\a\ 2 and a(3* + a* (3. 

In the case of the B — > D*~X C and -B — ► £)W~Jf u decays, similar considerations lead to 
the momentum spectra in the B rest frame 

dV iB^D;~X c ) = J^\a\ 2 lPD * f 



d \PD* s \ " '* '"' 2Km B E D *J\p D *\ 2 + m 2 c 



x 



f(C + c ) (2m B E 2 Dt C + c + m B m 2 Dt C + c - 3m 2 D ,E D *) - ($f - Cj] , (24) 



r^-^r s G\ [Tr T 2 2 2 \p Ds 



-(b - d;x u ) = -^\ Vub v; s \ 2 aif 2 D 



d\Pn s \" s 4W° csl UUa E Da J\ VDa \* + ml 
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X 



dr 



d\p D . 



where 



G 



f(C + ) (2m B E 2 Ds C + - m B m 2 Ds C + - m 2 D E D! 

2 



\v ub v; s \ 2 ajf D 



Aumi 



E D *J\v Dt \ 2 + ml 



x 



/(£») (2m B E 2 Dt C + u + m B m 2 m C + u - 3m 2 D ,E D; ) - (C + u - C 



(25) 



(26) 



E D * 3 ±^\PDj\ 2 + ml 
m B 



(27) 



For a longitudinally polarized D* we obtain 



d\p Dl 



-[B ^D*-(X = 0)X c 



f. 



IK 



a 







2 






2 + m 2 



X 



/(e; c ) (2m B \p D * \ 2 C + C + m B m 2 Dt C + c ~ m 2 D ,E Dt ) - (C + c - C 
For a transversely polarized _D*~ we have 
dT T 



(28) 



d|px>: 



-[B ^ D* s -(\ = ±l)X c 



7rm B 



\Q\ 



\PDi 



e d*,\ \Pd*\ 2 + m 2 c 



f(c + c ) (m B e + c -E D * - u*+ c ^c 



(29) 



From Eq. (2^) one can easily obtain the hadronic invariant mass spectrum for X c in 
B — > D~X C using the relation 

dT 'B^D 7 X c )= mx ; ED ° dT A B^D-X C ). (30) 



dm 



m B \p Ds \d\p Ds 



Similar formulas hold for the hadronic invariant mass spectra in B —>■ D* X c and B 

The total decay rates can be obtained by integrating Eqs. (|2"2"D, (|2~^)-(PB]), (^) and (| 
with respect to |p n (*)| over the entire kinematic range 



< |p D w| < J[{m B + m D (*)) 2 -m 2 x ][(m B - m w ) 2 - m\ ]/(2m B ), 



(31) 



where mx q . is the minimum value of the invariant mass of the hadronic system X q . 

Measurements of the ratio of the decay rate for the longitudinally polarized D*~ to 
the total decay rate for B — > D*~X C have a good potential for testing the theory, since 
the theoretical uncertainties partly cancel in the ratio and most of experimental systematic 
errors also cancel. 



7 



The leading effect of the initial B bound state is encoded in the distribution function 
/(C)- Since it is a nonperturbative QCD object, /(£) has not yet been completely determined. 
Nevertheless, several important properties of it are known in QCD P^JI7| . The distribution 
function is exactly normalized to unity, Jq d£,f(£) = 1, because of the conservation of the 
6-quark vector current by the strong interactions. For a free 6-quark, b(y) = e~ WPb b(0). 
From Eq. (p0|) it follows then that in the free quark limit 

MC = SU-— )■ (32) 



m B , 

The mean <£> = jj dtfftt) and the variance a 2 = ((£ - (£» 2 > = ft - (£» 2 /(£) 
characterize the location of the "center of mass" of the distribution function and the square 
of its width, respectively. They specify the gross shape of the distribution function. They 



have been estimated using the heavy quark effective theory. The results are ||16| , |17 

2 -(A x + 3A 2 



<0 = ^ 
m B 



Qmi 



a 



2 A; 

3m 2 , 



(33) 
(34) 



1 B 

where Ai and A 2 are the HQET parameters, which are defined as 

Ai = -^—(B\h v (iD) 2 h v \B), (35) 



2m b 
1 

12ms 



A 2 = -^—(B\h v g s G^h v \B). (36) 



Both of them are expected to be of order Aq CD . These imply that the distribution function 
is sharply peaked around m^/ms with a width of order AqcdA^b- This is in line with the 
expectation that the b quark should carry most of the momentum of the B meson. 

Another important property of the distribution function, discussed previously, is univer- 
sality. This implies that the distribution function can be extracted from measurements of one 
process, and then applied to predict others in a model-independent way. Indeed, the shape 
of the distribution function can be directly extracted from B — > tv(X q plj or B — > ^yX s fl7fl 
decays. This could eventually improve the precision of theoretical predictions. 

For numerical analyses we take the following parameterization for the distribution func- 
tion Eg 

/«)= v- ( y+v (37) 

where N is a normalization constant which guarantees Jq 1 d£f(£) = 1. The four shape 
parameters, (a,b,c,d), obey all the known constraints on the distribution function. Once 
the parameters c and d are given, the parameters a and b can be fixed by comparing with 
(£) and a 2 given in Eqs. (^) and Unfortunately, we do not know the values for c and 
d at present. We shall take c and d to be free parameters and vary them to gain an idea of 
how the momentum spectra for and branching fractions depend on the shape of /(£). 
Eventually, the shape of the distribution function can be determined from experimental 
data. The parameterization (|3"T|) facilitates the comparison to data. 
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IV. FREE QUARK LIMIT 



In order to understand the bound state dynamics of the B meson in the semi-inclusive 
B meson decays, we should consider the decay rates for the free quark decays, b — > D^~q, 
which give us a reference point. 

In the free quark limit, the distribution function becomes a ^-function. Substituting 
the free 6-quark distribution function Eq. fl32[) in Eqs. (p2|) [or, most easily, Eq. (|2l|) 1 and 
(f24|)-(|26|), we reproduce the corresponding total decay rates for the free quark decays in the 
rest frame of the h quark: 

m lf -™ 2 D s ( m b + m l)] + \P\ 2 { m l + m l ~ m D s ) 
m 2 c +m 2 D J}, (38) 



nb^D;c) = ^\p Ds \{\a\ 2 [(ml- 

+ (aj3* + a* f3)m c {m 2 — 



Y(b - D*-c) = J^\ a \ 2 \p Dt \[(m 2 b - m\) 2 + m 2 m {m 2 + m\ - 2m 2 m )}, (39) 
r(6 - Dju) = ^jV ub V;fa 2 J 2 Ds \p Ds \[(rn 2 b - m 2 u f - m 2 Ds (m 2 b + m 2 u )] } (40) 
r(6 - D*-u) = ^\V^V:fa 2 J 2 Dt \v D *J{m 2 - ml) 2 + m 2 D ,{m 2 b + m 2 u - 2m 2 D ,)]. (41) 



Similarly, substituting Eq. (^) in Eqs. fl28| ) and (p9|), we reproduce the correspond- 
ing decay rates for b — > D*~c with the longitudinal and transverse polarizations of D*~ , 
respectively, 



T L [b -> D*~(\ = 0)c] = J^\ a \ 2 \p D ,\[(mi - m 2 c ) 2 - m 2 Dl {m 2 b + m 2 c )}, (42) 



D t l„|2i I r/' 2 2\2 2 / 2 i 2\ 



T T [6 - D S *-(A = ±l)c] = ^%|a| 2 | Pl)| |K + - mi,.). (43) 

It follows from Eqs. (|3^) and fl4~2"D that the ratio of the decay rate for the longitudinally 
polarized D*~ to the total decay rate for b — > D*~c 

T L _ 2m 2 Dt (ml + m 2 — m 2 D ,) 

^ D * s C) = 1 ~~ (m 2 - m 2 c ) 2 + m 2 D% (m\ + m 2 c - 2m^ ) ' (44) 

Our results for the part of the free quark decay rate involving tree diagrams agree with 
Refs. |],|2|]. It should be stressed that the reproduction of the free quark decay results by 
taking the free quark limit of the leading twist results shows consistency of the light-cone 
expansion. 

All the momentum spectra for in the free quark decays b — > D^~q are a discrete 

line localized at 

IPdH = \J[(m b + m D M) 2 -m^mb-m^)) 2 -m 2 ]/(2m b ) (45) 
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because the momentum of in the two-body decay is fixed kinematically. 

As shown in the last section, at the hadron level the momentum spectrum for D^*^ in 
B — > D^~X q spreads over the kinematic range given in Eq. (|3TD , due to the fact that the 
hadronic invariant mass is changeable. However, since the heavy b quark in the B 
meson is nearly free, most of the spectrum remain around the free quark decay location of 
([451) ■ This is the phase space region where light-cone dominance occurs. 

For measurements of B — > D)*^~X U one can apply the cut on the 
momentum above the limit for charm production in the decay, |p n w| > 

^[{rriB + m D (*)) 2 ~ m D}[( m B — m D { * ) ^ 2 ~ m £>]/ t° suppress the B — > D^~X C back- 
ground. Since this cut is well below the free quark decay location of ((4"5|), the majority of 
the momentum spectrum in B — > Dy>~X u pass the cut, as we shall see. 



V. NUMERICAL ANALYSES 

For numerical analyses, we must specify the values of the input parameters. We use 
G F = 1.16639 x 1(T 5 GeV~ 2 , m B = 5.279 GeV, m Da = 1.969 GeV, m D * = 2.112 GeV, 
m D = 1.87 GeV, and m n = 0. We neglect the B + , B° lifetime difference and use r B = 
1.6 x 1CT 12 sec. For quark masses, we use m b = 4.9 GeV, m c = 1.5 GeV, m s = 120 MeV, 
and m u = 0. As for the CKM matrix elements, we take \V u b\ = 0.0035, |V CS | = 0.9742, 
\V c b\ = \V ts \ = 0.04, and \V tb \ = 0.9992. Since the phase u ud = &rg(-V* s V t * b V cb V ts ) is very 
small, we assume coscow = 1. For the decay constants, we take fz> s = 264 MeV from 
the measurements of the leptonic branching fractions Br(D s — > tv) and Br(D s — > fiu) f23fl 
and assume fn* = fn. We take mv = m c for B — > D{*^~X C and m x „ = for 

"8 s c min * u min 

B — > D(*'~X U . Finally, we consider two very different shapes of the distribution function 



parameterized in Eq. (37), corresponding to two sets of parameters: (i) preset c = d — 1, in 
that case a = 0.9548 and b = 0.005444, which are inferred from the known mean value and 
variance of the distribution function given by Eqs. (Q) and (|34[) using the HQET parameters 
Ai = —0.5 GeV 2 and A 2 = 0.12 GeV 2 ; (ii) preset c = d = 2, in that case a = 0.9864 and 
b = 0.02557 inferred from the same mean value and variance of the distribution function 
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We first compute the branching fractions for B — > D^~X q . It turns out that the 
branching fractions are not sensitive to the detailed shape of the distribution function. The 
difference between the branching fractions obtained using the parameter sets (i) and (ii) 
for the distribution function is negligible. The results for the direct production of in 
semi-inclusive B decays are 

B(B -> D~X C ) = 4.0%, 
B(B -> D*-X c ) = 5.4%, 
B(B -> D;X U ) = 4.8 x 10~ 4 , 

B(B -> D*-X u ) = 6.2 x 10" 4 . (46) 

The calculated branching fraction for B — > D*~X C is in agreement with the measured one by 
BABAR quoted in the Introduction. Because D*~ decays to D~ r y or D~7i°, we can compare 
the sum of both channels, B(B -> DjX c ) + B(B -> D*'X C ) = 9.4%, from Eq. (|BJ) with the 
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measured branching fraction for B — > DfX quoted in the Introduction. The comparison 
shows an encouraging agreement. 

We find that the contribution of the penguin amplitude decreases the branching fraction 
for B — > D~X C by 3% and the branching fraction for B — > D*~X C by 6%. So the penguin 
contribution is small but non- negligible, especially for B — ► D* _ X C . 

Comparing with the calculations in the free quark decays b — > D^~q, we find that all 
the branching fractions for B — > D^~X q are enhanced by the initial bound state effect by 
about 5%. 

The ratio of the decay rate for the longitudinally polarized D*~ to the total decay rate 
for B — ► .D* - X c is also insensitive to the detailed shape of the distribution function. The 
bound state effect increases the ratio by 2% compared to the free quark decay. We obtain 

^(B - D* s -X c ) = 0.66. (47) 

Now we turn to the differential decay rates. Figures 1-4 show the momentum spectra 
for £>W- i n B -> We choose to use two very different shapes of the distribution 

function to calculate the spectra. In contrast to the integrated decay rates, all the spectra 
depend strongly on the shape of the distribution function. Hence measurements of the 
spectra can provide useful information about the nonperturbative distribution function. 

As Figs. 3 and 4 make clear, the majority of the momentum spectrum for D^*'~ in the 
decay mode B — > Dy'~X u lie above the kinematic limit for the dominated decay mode 
B -> DW-X e , which is |p D J > 1.81 GeV for B -> D~X U and |p D *| > 1.73 GeV for 
B — > D*~X U . The fraction of spectrum above the charmed limit is not sensitive to the 
shape of the distribution function. We find that about 98% of the spectrum goes beyond 
the charmed limit. Thus, the kinematic cut on the momentum of provides an efficient 
discrimination between B — > D(*>~X U decays and B — > D(*'~X C decays. 



VI. DISCUSSION AND CONCLUSION 

We have calculated the total and differential B — > D^~X q decay rates assuming fac- 
torization. We have treated the initial bound state effect using the light-cone expansion 
and the heavy quark effective theory, generating a theoretical description of this effect from 
QCD rather than the phenomenological model 0. The resulting branching fractions for 
are found to be in agreement with the measurements. This suggests that 
factorization is a fair approximation for these decay modes. However, there are still consid- 
erable uncertainties in the theoretical calculations and experimental measurements. Theo- 
retical uncertainties arise from the input parameters, the 6-quark distribution function, the 
higher twist corrections, and the radiative corrections. The uncertainty in the decay rate 



from the error on the decay constant fo s alone |23] is already quite large (~ 30%). 

For the determination of \V u b\ from B — > D^~X U , additional theoretical uncertainty 
results from the underlying factorization hypothesis itself. It is hard to quantify this uncer- 
tainty. Therefore, \V u t\ determined from the semi-inclusive nonleptonic decay B — > D(*>~X U 
will not be competitive with that determined from the inclusive semileptonic decay B — > 



viX u [24]. Nevertheless, B — > X u will provide a useful independent measurement of 
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1 14b I with very different systematic errors. On the other hand, using \V u f,\ determined from 
the inclusive semileptonic B decay, the range of validity of calculational tools for addressing 
QCD can be tested and established in B — > D^~X U , 
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FIG. 1. Momentum spectrum for D~ in B — > D~X C . In Figs. 1-4, the solid curves are obtained 
using the parameter set (i) for the distribution function; the dashed curves are obtained using the 
parameter set (ii) for the distribution function. 
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FIG. 2. Momentum spectrum for D*~ in B — > D*~X C . 
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FIG. 3. Momentum spectrum for D~ in B — > D~X, 
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FIG. 4. Momentum spectrum for D* in B — > D* X„. 
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